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Curve Fit Method, 50% Transmitted Light ---------- 50 
Two Angle Method, 50% Transmitted Light ---------- 51 
Curve Fit Method, 30% Transmitted Light ---------- 52 
ive mete ethed, S0s)ircansmitced Light —==—-—--—--—- 23 
2% Aluminum Propellant, Less Than 90% 

Attenuation ----- ee nn a 54 
HNeie)omteiZocmerooe |Panteml2 OCE., Exhaust ------- 252) 
ose reine bopellant ys 1S Oct., Exhaust -—------ 56 
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4.8% Aluminum Propellant, 18 Oct., Exhaust ------- SF 
4.8% Aluminum Propellant, 23 Oct., Exhaust ------- 58 
4.8% Aluminum Propellant, 23 Oct., Exhaust ------- 59 
2% Aluminum Propellant, 30 Oct., Exhaust --------- 60 
2% Aluminum Propellant, 30 Oct., Exhaust --------- 6 
4.8% Aluminum Propellant, 19 Nov., Exhaust ------- 62 
4.8% Aluminum Propellant, 19 Nov., Exhaust ------- 63 
2% Aluminum Propellant, 25 Nov., Exhaust --------- 64 
2% Aluminum Propellant, 25 Nov., Exhaust --------- 65 
2% Aluminum Propellant, 26 Nov., Exhaust --------- 66 
2% Aluminum PEOseiiane 26 Nov., EXRAUSt #2 67 
SEM Eval., 2% Al Propellant see 68 
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The importance of aluminum oxide (AL 503) particle size 
produced in solid-propellant rocket motors has been recog- 
nized since aluminum was first introduced as a fuel additive. 
Aluminum added to a propellant will increase the thrust 
produced and aid in suppressing high-frequency combustion 
instability. Although the performance improves with aluminum 
addition, two phase flow losses adversely affect specific 
impulse efficiency. Specific impulse efficiency is the product 
of several component efficiencies, i1.e., chemical-kinetic, 
combustion, two-phase flow, two-dimensional flow, submergence 
and throat erosion. In aluminized propellants, the two-phase 
flow loss is one of the more significant losses. Two-pnase 
flow losses result from particle-gas thermal and velocity lags. 
The latter strongly depends on particle size. Parametric 
studies involving small motors, with throat diameters less 
than one inch, at AFRPL [Ref. 1] and CSD [Ref. 2] have 
revealed that a 1 micron change in particle diameter results 
in a change of 1% in Isp. The effective use of particle 
damping for suppression of combustion pressure oscillations 
is also strongly dependent upon the particle size distribution. 

One of several analytical computer programs used to pre- 
dict performance of a solid-propellant rocket motor is the 


SPP (Solid Performance Program) [Refs. 3,4]. As the above 
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discussion indicates, it is of critical importance to accurately 
predict particle size. This is especially true for particle 
damping with minimum smoke propellants which typically con- 
tain less than 2% metal by weight. The 1975 version of the 

SPP [Ref. 3] contains a particle sizing empirical model 
developed by Cohen. This model was based on correlating 
collected exhaust samples. After a critical review, this model 
was considered inadequate. A new model was therefore required 
for a new SPP version, titled Improved SPP [Ref. 4]. Inorder 
to meet SPP accuracy goals, this model had to predict particle 
size within +10%. Unfortunately, there are no adequate 
theoretical models relating particle size to propellant and 
motor parameters. This necessitates the use of empirical 
methods for correlating particle size. 

There are two methods that can be employed to develop an 
empirical relation. One is to correlate particle size data to 
propellant and motor variables, as was done by Cohen. The 
other method is to calculate particle size based on a critical. 
Weber number for a maximum stable drop size [Ref. 5]. After 
considering the relative merits of each method, Hermsen [Ref. 
5] devised an empirical model based on collected particle size 
data for the improved SPP. This model was correlated to a much 
broader data base than was previously available to Cohen. 

The particle diameter correlated was Daas the mass-weighted 


diameter, given by: 


iuak 


_ 4 3 
Dgs 7 
Smith [Ref. 6] has shown that D3 was the most appropriate 
average to use in predicting two-phase flow losses. 


Hermsen's empirical relation [Ref. 5] is presented 


below: 
oe = 3.6404 p22932 (4 ~@7 (9008163EPT) 
ie 
where: 
Dae maSss-weighted average diameter, Ha 
D, = nozzle throat diameter, in. 
E = AL,03 concentration in chamber, g-mol/100 g 
P = chamber presSure, pSia 
ze = average residence time, msec 


The predicted Da3 iS now applied to the expression below in 


order to calculate the two phase flow loss [Ref. 7]. 


& 6: 


- _ 5 (eee see 6 

ETATP = 3 ee /P E D. ) 

where: 

ETATP = impulse loss effect, % 

S = mole fraction of condensed species, moles 
of condensed species per 100 grams of 
mixture 

D., = Dy3 TASMPeroms 


eZ 


P = average motor chamber pressure, pSia 


E = nozzle expansion ratio at ignition 
COV ited ons 

D, = nozzle throat diameter at ignitiion, 
inches 

= = Naperian base 2.71828 

C( ) = coefficients that are a function of Dd, 
and De 


The two-phase flow loss is applied with the other aforementioned 
efficiencies to a theoretical Isp. The Isp loss from boundary 
layer effects is then subtracted from the corrected theoreti- 
cal Isp to obtain delivered Isp. 

Collecting exhaust produces: is feasible only for small 
motors and the techniques used often result in variation in 
the measured size. Validation of the SPP model must be based 
on data collected within the motor environment. Light 
scattering measurements are well-suited for this application. 
These measurements are more easily made if the particles are 
large compared with the wave length of the transmitted light 
and the transmittance is greater than 90%. The former allows 
the accurate application of Fraunhofer diffraction, whereas the 
latter specification is necessary in order to satisfy single 
scattering requirements. 

The Naval Postgraduate School has conducted a series of 
investigations to determine the applicability of light scatter- 


ing measurements to the solid rocket motor exhaust products 


iS 


[Refs. 8,9,10,11]. This investigation was directed at modi- 
fying the apparatus and techniques developed in the previous 
efforts in order to more accurately determine particle size. 
Validation of the apparatus was also a goal. The SPP particle 
prediction model can then be verified by a comparison of the 


predicted size to the experimental data. 


14 


If. THEORETICAL BACKGROUND 


to 


The most general theory of light scattering was developed 
by Mie and is presented by Van de Hulst [Ref. 12]. Mie's 
theory, although not restrictive in its assumptions in regard 
to particle size and refractive index, is however very com- 
plex to utilize for data reduction. Dobbins, et al. [Ref. 13] 
and Powell, et al. [Ref. 14] derived similar expressions for 
intensity ratios at two forward scattering angles based upon 
an Upper Limit Distribution Function (ULDF). and Fraunhofer 
diffraction. The ULDF was initially proposed by Mugele and 
Evans in Reference 15. Although the expressions developed by 
Dobbins and Powell are extremely useful in relating intensity 
ratios to Door the volume-surface mean diameter, these models 
are still somewhat cumbersome to use for data reduction 
purposes. 

Particle sizes found in solid-propellant rocket motors 
are often bi-modal. Particles in the larger mode are the most 
important for two-phase flow loss calculations and are gener- 
ally much greater in diameter than the wavelength of inci- 
dent light. Scattering by these large particles is described 
by Fraunhofer diffraction. Hodkinson [Ref. 16] and others 
{[Refs. 13,14] have demonstrated that diffraction theory 
(within the forward lobe) can be used to approximate the Mie 


theory. Buchele in Reference 17 summarized the experimental 


iS 


techniques for determining the particle size by measurements 
of diffractively scattered light. Buchele presents a function 
which closely approximates the curves given in References 13 


and 14. 


1(8) = EXP[-(.5708)*] 


where: 


I(§) = ratio of intensity of scattered light at 
some angle (98) to the intensity of 
scattered light at theta equal to zero 
degrees 

a, = m7D/A is the particle size parameter for 
diameter D2 and wavelength of light 
lambda (Aj); 


The above expression was used in the current study for data 
reduction to obtain D25- Unfortunately, D3 is the mean 
diameter used in the SPP. Therefore, D5 must be related to 
Dy 3 either by analyzing exhuast samples or by approximate 
analytical means. 

A final important point is that the above expression is 
valid for a < 3 (the center lobe) 9) ror vthete pected eran 


of theta and D this proves to be a very mild restriction. 
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Tit. BXPERIMENTAL APPARATUS 


A. ROCKET MOTOR 
The rocket motor used in this investigation was the same 
as the long motor described in Reference 11. The motor com- 


ponents are presented in a photograph in Figure 3.1 and the 


installed motor 1s shown in Figure 3.2. A schematic is also 
presented in Figure 3.3. The propellant grain dimensions re- 
mained the same as in the previous effort [Ref. 11]. The 


grain was two inches in diameter, one inch in length and had 
a web thickness of .725 inches. The grain was cylindrically 
perforated and allowed to burn also on the aft end. The 
igniter was composed of BKNO,. Ignitiion of the igniter was 
achieved by heating a resistance wire with a 12 VDC source. 
Some of the exhaust products were collected in an eight 
inch diameter stainless steel exhaust tube. The products 
were collected from both the side walls of the tube and the 
end cap. Another particle sample was collected from the walls 
of the combustor, near the nozzle entrance. Both samples were 
then cleaned and mounted for observation using a scanning 


electron microscope. 


B. LIGHT SCATTERING APPARATUS 
The apparatus shown in Figures 3.4 and 3.5 was basically 


the same as used in Reference 1l with some exceptions. The 


ie 


light source used for the nozzle exhaust was a 8-miliwatt, 
helium-neon laser. The nozzle entrance also eed a 8-milliwatt 
helium-neon laser. The lasers were mounted on two parallel 
optics benches, such that one beam could pass through the ex- 
haust and the other beam could pass through the motor cavity 
adjacent to the nozzle entrance. A beam expander/collimator 
was used with the Nelium-neon laser in the exhaust. Each beam 
was intercepted by a stop placed directly in front of the 
narrow pass filter. The narrow pass filter was used to eliminate 
extraneous light. The stop and narrow pass filter were placed 
30.5 centimeters from the motor centerline for the exhaust 
beam and 11.5 centimeters from the motor centerline formmed= 
motor cavity beam. The exhaust beam narrow pass filter had 

a Giameter of 5.08 centimeters and the motor cavity narrow 

pass filter had a diameter of 2.3 centimeters. Both beams 
passed through the centers of these filters. In the present 
arrangement of the apparatus this resulted in a maximum angle 
for collecting scattered light of 4.8 degrees in the exhaust 
and 6.0 degrees in the motor cavity. The placement of the 
optics closer to the motor centerline allowed for a greater 
maximum angle to be recorded than in the original design. 
Directly following the narrow pass filters were condensing 
lenses. Both paths used a lens with a diameter of 5.08 
centimeters. These lenses had a focal length of 50 centimeters 
and were utilized to focus scattered light onto a linear array 


of photodiodes. The array consisted of 1024 silicon photodiodes 
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with .025 millimeter spacing. The apparatus in Reference 1l 
positioned the lenses such that only angles from 0 to .04 
mitans COUld be observed. In actuality, however (as a conse- 
quence of diffraction around the beam stop and beam spreading), 
the range of useful angles where only scattered light occurred 
was from .02 radians to .04 radians. In the current arrange- 
ment the diode array was placed 1.0 centimeters below the main 
beam and the collecting tube diameter was enlarged. This 
allowed a range of angles from .02 to .07 radians to be 


observed. 


C. DATA ACQUISITION AND DATA REDUCTION 

The data acquisition system 1s presented in detail in 
Reference 10. The system controller was an HP 9836S computer 
and the HP 6942A multiprogrammer vrovided A/D conversion and 
Storage. The data acquisition system was programmed to take 
data at a specified point in the motor firing by inputting a 
trigger pressure and a time delay after the trigger pressure 
was reached. The number of diode scans for the exhaust anc 
motor cavity were eight and four respectively (limited by 
the current data acquisition system memory size). 

Two methods were used in References 10 and 11 for data 
reduction, one of which was altered for this study. The modi- 
fied method pertained to the iterative graphics technique. 

In References 10 and 11 the experimental data were normalized 
and compared to a theoretical profile that was also normalized 


to some assumed scattered light intensity at zero degrees 


ile, 


for a specified D A "best-fit" of tne experimental profile 
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to the theoretical one yielded D The current investigation 


Be 
determined: that this method resulted in ambiguity in deter- 
mining D5 aS a consequence of assuming a centerline intensity 
for the scattering profile. The method used in the present 
study was one in which the experimental data were not altered. 

A theoretical curve was generated for any specified value of 

D35 based upon a selected initial angle (05) and the correspond- 


ing intensity at that angle (I,) from the unaltered experimental 


al 
data. The approximate intensity profile for a polydispersion 


given by Reference 17 is: 


5 2 2 2 
ae I, ExXP-{(6, Snow nero? Diy ee} 


Specification of a series of values for 0. results in 


values for I, and, therefore, an I vs 6 theoretical profile 


2 
that passes through the experimental data at Tis ea abit 

this curve does not fit the experimental data, then other values 
OF D35 can be chosen and a new theoretical curve generated. 

A "best-fit" method is subsequently used to determine D35- 
In summary, the method used for obtaining the mean particle 
Size 1s based upon fitting an approximate (although quite 
accurate within specific angle restrictions) theoretical 
profile for a polydispersion to the unaltered profile obtained 


experimentally over a small range of scattering angles. The 


modified data reduction program is presented in the Appendix. 
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The second method was similar to the first and also used 


the above expression solved for D325: 


Zee 2 
D = ~Ln(T,/I,) (A/7.577) X68, - 8) 
The procedure used was to first select a starting angle, O57 


2 


from the experimental data. This was repeated with increasing 


along with a 65/84 ratio. A value of I, at 8, was obtained 


values of 6, and then again with different values of 64/8). 


Z 
The results were plotted graphically as particle size vs 
scattering angle for each angle ratio. A nearly constant 


value of D3. vs 8 indicated a good correlation. 
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Description 


Outside Diameter 


(dew iieh)) 


Length (ln, inch) 


Throat Diameter (dt, inch) 


Slope (converging) 
Angles yy, Gdeguee) 


TABLE I 
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Figure 3.1. Motor Components 
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LEGHT SOURCE SPEC#EHEATION 


Helium-Neon Laser (Exhaust) 


1. Manufacturer Spectra-Physics 
2. Model: a7 

3. Type: He-Ne Class IIIB 
4. Output Pewer: 8 mWatt 

5. Beam Diameter: «92 in 

6. Beam Divergence: -o/ fmaad. 


Helium-Neon Laser (Motor Cavity) 


I; Manufacturer. Uni-Phase 

2. Model: 3035 

3. Type: He-Ne Class IIIB 
4, Output Power: 5 mWatt 

5. Beam Diameter: ; 61) aim 

6. Beam Divergence: 1. 0DURimmaae 
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fee = LNTRODUCTION 

The purpose of this study was to improve and validate 
the accuracy of the apparatus used for particle sizing in 
the rocket motor environment. Also, the experimental data 
will provide a basis to verify the particle size prediction 
mage] Of the SPP. 

The apparatus had been calibrated by Harris [Ref. 10} 
and Kertadidjaja [Ref. 11] using glass and polystyrene spheres 
in a suspension of water. Kertadidjaja [Ref. 11] also 
attempted to demonstrate the applicability of the apparatus 
to the combustion environment. This investigation continued 
the earlier studies with several refinements. Three propel- 
lant types were used in this investigation. All three utilized 
a GAP binder with ammonium perchorate (AP) oxidizer. Each had 
varied aluminum content, with the percentages present being 
03, 2%, and 4.8%. The size of the aluminum powder was 20 
microns. The major ingredient compositions for the 4.8% 
aluminized propellant were 14.7% GAP, 45.7% 200 um AP, 24.6% 
fereum Ae, 6.5% TEGIN and 4.8% 20 um AX . The 0% and 23% 


aluminized propellants kept the same solids loading. 


B. SYSTEM CALIBRATION 
As previously mentioned, the apparatus was calibrated by 


both Harris [Ref. 10] and Kertadidjaja [Ref. ll]. The 


Ze 


apparatus was re-calibrated in this study as a result of the 
modifications that were made. The standard calibration proce- 
dure described in detail by Hare [Ref. 10] remained the 
same. The laser beam passed through a glass container holding 
a suspension of polystyrene spheres and water and was initially 
positioned on the number one diode. The linear diode array 
was then moved vertically down, with the aid of a dial indica- 
tor, such that the number one diode was one centimeter below 
the main beam. Since the diode array was 2.56 centimeters 
long and was 50 centimeters from the focusing lens, this 
arrangement allowed for a range of meaSuring angles from .02 
to .07 radians. Therefore, the scattered light striking the 
diodes would not be affected by the main transmitted beam or 
the light diffracted by the beam stop as in the previous two 
studies. This arrangement was used both for the exhaust beam 
and the motor cavity beam. 

Two polydispersion samples of polystyrene spheres, with 
D5'S of 4.8 microns and 10.2 microns, were used for the cali- 
bration. These samples were chosen because the D.5 was known 
and the expected D5 of the actual aluminum oxide particles 
during a motor firing was between 1 and 15 microns. The 
results are presented in Figures 4.1, 4.2, 4.3 and 4.4. These 
results were encouraging in regards to validating the modifi- 
cations and the method for data reduction. Calibrations were 
also conducted within the motor cavity with the laser beam 


and scattered light passing through the motor windows. Figures 
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4.1 and 4.2 are representative of the results obtained. This 
indicated that the motor windows did not adversely affect 
scattering profile and, therefore, ee the results. 

An initial point of concern was the narrow range of measur- 
ing angles which were employed in the present apparatus. 
Since large particles diffract most of the incident light at 
small angles and very small particles diffract most of the 
incident light at large angles, the range and minimum angle 
employed could not have resulted in sufficient sensitivity to 
particle size for the expected range of particle sizes. Also 
shown in Figure 4.1 are two additional theoretical profiles 
mer D of 10 and 8.3 microns. It was apparent from the sensi- 


e2 


tivity of the slope change with D that the mean particle 
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size could be determined with an uncertainty less than +.5 
mubcrons . 

The "upsweeps" in the curves at larger angles using the 
two angle method were the result of the increase in slope of 
the scattering profile due to tne effect of truncation by 
the apparatus. These parts of the curves remain only to 
provide a reference to the angle ratios. 

Another informative point .in apparatus design is given 
by Hirleman [Ref. 18] in regards to placement of the focusing 
lens. Hirleman notes that truncation of the larger scatter- 
ing anbles, which is a function of the placement of focusing 


lens, will cause a biasing against the small particles. The 


expression given in this discussion is: 


one 


Z = (md/7.66A)D 


Cc 
where: 

Zo = critical distance in front of the focusing 
lens where significant energy diffracted 
by particles less than diameter D is 
truncated by the focusing lens 

IN = wavelength of the incident light 


For a given Zoe all particles that are.greater than D in size 
will have 84 percent of their diffracted energy pass through 
the focusing lens. In the exhaust configuration of the present 
apparatus Z was equal to 30.5 centimeters, d was 5.08 centi- 
meters and A was .6328 microns. This yields a D of 9 microns. 
Thus, according to this criterion, some of the scattered light 
from particles below 9 microns is truncated and, therefore, 
results ina bias towards larger particles. This type of 
biasing was small in the present apparatus down to particle 
sizes of approximately 3 microns (the smallest particles in 
the 4.8 mine won D35 polydispersion) as evidenced by the good 
calibration results is shown in Figure 4.3. However, the 
apparatus could be improved by either moving the lens 
closer to the exhaust (at the risk of apparatus damage from 
the exhaust jet) or by using a larger diameter lens. 
Transmittance effects were also investigated in conjunc- 
tion with the calibration. As previously mentioned, trans- 
mittances above 90% are required to meet single scattering 


requirements. Calibrations were conducted using polystyrene 


Sie 


spheres with a D5 of 10.2 microns. Scattering data were 
recorded for transmittances of 85%, 70%, 60%, 50%, and 308%. 
meeresulte are presented in Figures 4.5 through 4.14, and 
summarized in Table III. Both the curve fit and the two- 
angle methods resulted in very good agreement with the actual 
mean particle size with transmittances as low as 50%. Evena 
transmittance of 30% did not severely affect the measurement. 
These results indicate that particle sizing methods based upon 
Single scattering theory can effectively be used in a multiple 
scattering environment with transmittances significantly less 
imam 9OR. 

A test was then conducted to determine the transmittance 
in the actual motor exhaust uSing the 2% aluminized propellant. 
The procedure was the same as the one used in Reference ll. 

The transmitted beam was directed onto the diode array (with 
the appropriate filters to reduce the intensity) and the beam 
profile was recorded. The motor was then fired and the profile 
waS again recorded. The results are presented in Figure 4.15. 
The results indicated a transmittance greater than 90%, with 
some beam spreading. There was also no apparent beam shift. 
Therefore, transmittance leoseevascsdecrermined not to be a 
debilitating factor in the motor exhaust jet. 

The calibration curves have an apparent limit of approxi- 
mately .04 radians for which scattered light can be recorded 
(for example, see Figure 4.3). However, it was stated that 


the apparatus was designed to record angles up to .07 radians. 
| 
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This is due to the water-glass interface having an effective 
index of refractionwet 12729. Thus ede scattering angles are 
corrected for this high index of refraction, resulting une 


corrected maximum angle of approximately .04 radians. 


C. NOZZLE EXHAUST MEASUREMENTS 

One test was conducted with a nonmetallized propellant. 
The results, presented in Figure 4.16, show no Sscattermace 
Therefore, the scattering observed by Kertadidjaja [Ref. 11] 
at smaller angles for a similar test was probably the result 
of beam broadening and diffraction around thembeeam store 

Three tests were conducted with both the 4.8% and 2% 
aluminized propellants. The motor firing data and a summary 
of the particle size measurements are presented, respectively, 
in Tables IV and V. From light scattering measurements each 
propellant showed a similar result, one test with a D5 Om 
approximately 5 microns and one with a larger value of D3. 
(7-12 microns). This test-to-test variation probably resulted 
from the very short burn times and the corresponding lack of 
a well-defined steady-state condition. These results are 
presented in Figures 4.17 through 4.28. Typical S.E.M. photo- 
graphs are shown in Figure 4.29. The light scattering results 
were quite consistent with those obtained from the S.E.M. 
evaluation. Although the sample sizes from which the S.E.M. 
evaluations were based were relatively small (see Table VI), 


it still provided a good indication as to the validity of 
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the light scattering measurements. It should be noted that 


D351 and to an even greater extent D is strongly dominated 


Ae 
by the larger particles. One or two large particles can 
Significantly change D35- Future S.E.M. evaluations should 
utilize much larger sample sizes. 

The SPP equation for Da3 wilehn uS primarily a £Lunction of 
throat diameter) predicted much smaller values of Dy3 than 
observed in this limited set of data. The SPP model, however, 


is based on particulate data from motors with nozzle diameters 


greater than 1 inch. 


D. NOZZLE ENTRANCE MEASUREMENTS 

Light scattering measurements through the motor at the 
nozzle entrance were unsuccessful as a result of combustion 
light existing at the same frequency as the laser light. fThis 
was the case for both the helium-neon laser (A = .6328 um) 
and the argon laser (A = .488 hum). 

The analyses of the S.E.M. data for the nozzle entrance 
indicated a much larger D345 and Da3 as compared to the nozzle 
exit measurements. Therefore, a significant reduction in size 


fee abowto ave Occurred to the larger particles across the 


exhaust nozzle. 
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TABLE user 


EFFECT OF TRANSMITTANCE 
ON MEASURED PARTICLE SIZE 
(D35 from Manufacturers Data Equals 10.2) 


Transmittance Measured D 
(3) (microns) 
85 9.8 
70 76 
60 9.5 
50 9.4 
30 922 
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TABLE IV 


SUMMARY OF MOTOR FIRINGS 


Date of Wt of Press Press Burn Residence Dy 
Test Alum Po Max P Time Time (ar ) 
(3) (psig) (psig (sec) (msec) 
wee 31, 85 0 350 366 Beis) - 28 
eet 22, 85 4.8 Baz 328 ae Or ~28 
@et 23, 85 4.8 323 B59 Pigs ee. 3 aS 
See 18, 85 4.8 B26 340 208 lS a2 3 
Mew 19, 85 4.8 450 466 oS 14.0 2D 
wer 30, 85 230 296 2916 28 10.4 226 
Ney 26, 85 2d 408 467 20 14.4 . 28 
Ney 25, 85 20 566 57S 1 iG 27 B25 
Note: Pressure colume labeled P_ indicates pressure 


at which data was taken. 


a7 


Wt of 


Alum 


(%) 


4.8 


4.8 


4.8 


ZnO 


ZaO 


220 


Press 


(psig) 


Note; 


328 


333 


450 


296 


408 


566 


SUMMARY OF EXPERIMENTAL RESULTS 
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TABLE V 


D39 
S.E.M. 
Ext. /Mtr. 


52 6/LO.2 


136713 .9 


227 1454 


-/- 


eco/lOse 


D4 3 
S.E.M. 
Ext. /Mtr. 


7.3/24.7 
-/- 

19.7/23.9 

18.6/19.6 
-/- 


6. O/slse2 


3 
SPP Model 


2a 


250 


Ze 


PASS: 


ZieD 


ZoS 


(-) indicates sample was not taken for this run. 
Diameters are in microns. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The results of this study have shown that the measure- 
ment of D35 in rocket motor exhausts can be made using quite 
Simple light scattering methods. Apparatus calibration 
results and S.E.M. evaluations of the collected motor exhaust 
products support this conclusion. 

The SPP equation model underestimated D3 for this inves- 
tigation. This indicated that the model predictions cannot 
be extrapolated to small motors with throat diameters less 
than one inch. The light scattering method, however, can 
be used to validate the SPP model for the larger motors upon 
which it was based. 

In order to improve the signal to noise ratio in regards 
to the scattered light intensity, it is recommended that the 
laser beam be moved closer to the exit plane of the nozzle, 
where the particle density iS greater. The laser beam for 
this study was 7.3 centimeters aft of the nozzle exit. To 
resolve the combustion light problem it may be necessary to use 


eo ee ecenes laser, suen asa Cyn that has a different wavelength 


2 
than present in the combustion process. Also, other binder 
systems should be tested to determine if they significantly 


alter the light spectrum. 
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APPENDIX 


PROGRAM LISTINGS 


VS AEE ESSE ESE SE SE HE SESE SESE EE TE EEE EE 


VOX KERRY ERR EEM ENS RHC HE KELES LELEX REEL KEY 8 
VEX ERK EE RK EAE EEE RK E PLOTS RAW DATA Tee POC tet ecee re ets 2: 
VERRY R RAK EER MY RENEE AVEBACES BX BR YEY HELE HEXEN 
eecitatt tet tte ttt te: FILTERS {ECBO OR AGE 


LexXEXEX RR MRMERELK RX DETERMINES MEAN DIAMETER xx xe HEX KEREX ONE EY X 
VHRR ERREX HREM HEKHES RY TNTERACTIVE GRAPHICS x¥ 4H HHH EH HE RHEL EELS 
Hauer ereecee AND THE TWO-ANCLE METHOD #4 X4KEX CHEEK FRREK EX F 
Pireteteertisteres tT: Robert Kelly Harris BA €#xHHRE KEL KK CHEN AEE 
VERE EY ERR E ERE EL ELE Revised, John S Rosa AG ee HHREX EHEEH XEREXEK X 
BBE GRGOGBE OBIE EEE GEL 
OPTION BASE | 

COM /Hrdqaus/ Av2(1024) 

COR A2(1924) RUFFER 

COM /Gauss/ Tt(1024),L 

COM /Max/ M7 M5, Xt, Yt, Xm, Ya, Xx 

COM /Readata/ B,P,H,Q3$(20) ,Q4$(20) 278120) ,Y1(8192) BUFFER, Y2(8192) BUFFE 


COM /Two/ Av1(1024) MMI F 

DIM Scans(8) ,X11924) 

POM /Plots/ PIS{20) ,P2st2n) P3sl20) ,P4sl2o) 
INTEGER Sraf(t:t2480) BUFFER 


210 Choose: PRINT CHRS(12) 


ord 
259 
ON" 
24) 
oo4 
oo 
70 


CB) 


079. 


3Ng 
319 
324 
330 
340 
Sef 
350 
379 
380 


Old=9 

PRINT * TO LOOK AT ANY DATA FILE THE PROGRAM NEEDS SOME STARTING INFORMATI 
PRINS ** 

Pee = TO ACCOUNT FOR THE CHANGE IN WAVELENGTH IN THE MEDIUM" 
PRINT "* 

PRINT * ENTER THE INDEX OF REFRACTION OF THE HEDIUM" 

Pe UH lies 

PRINT " WATER = 1,353" 

PRINT “ AIR = 1.0" 

PRINT ° ESTIMATE OF EXHAUST = 1.1" 

REEP 

INPUT " THIS VALUE ADJUSTS THE COMPARISON CURVE TO THE MEDIUM", MI 


PRINT CHR$(12) 
PRINT * TO ACCOUNT FOR REFRACTION OF LIGHT AT BOUNDARIES KETWEEN® 
PRINT "” 


PRINT " THE MEDIUM AND AIR YOU MUST ENTER THE INDEX OF REFRACTION " 
PRINT "" 


70 


a aM a OF THE COMBINATION OF THE MEDIUM AND ITS BOUNDARY" 

410 PRINT * THIS VALUE IS APPLIED DIRECTLY TO THE DATA" 

420 PRINT ** 

43) PRINT * AIR=1.0 PLEXIGLASS & WATER = 1.39" 

440 PRINT ** 

450 PRINT “ESTIMATE OF MOTOR CAVITY COMBUSTION PRODUCTS & WINDOW = 1,22" 

460 PRINT ** 

47) PRINT * ESTIMATE OF ROCKET EXHAUST = 1.1 OR 1.0° 

480 INPUT ¥ 

49) PRINT CHR$(12) 

500 = INF(IT * ENTER LASER WAVELENGTH (HeNe=,6328 Ar=.489)",L 
S10 INPUT "ENTER THE RESIRED PLOTTING INTERVAL OF DIONE APRAY (2,4,6)",H 

Son  F=500 lnm FOCUL LENGTH OF OBJECTIVE LFNS 
Set Xtt)=1 

522 Tt )=,92 

539 Diod=.925 Inn DIODE SPACING MAY RE .03 

549 FOR I=1 TO 1023 

550 X(L+1 )=T+1 ICREATE AN ARRAY OF DIOTE NUMBERS 
50 THCT#LSATNC CL eDi od) /(F HM) 4T1 1) ICREATE AN ARRAY OF THETA 

57 NEXT | 

580 PRINT CHRS({2) 

590 CALL Displavi (Old) I JUST PRINTS INITIAL 2FKAR<S ON CRT 
600 IF Old THEN CALL Review(Avl(#)) ITHESE TWO LINES ALLOY FOR 

610 IF Old THEN COTO Gauss IREVIEWING REDUCED DATA FILES 

62)  REEP 

630 INPUT * INPUT TWO FILENAHES NOW (NO-PART PART)" O38, 04 

640 = INPYT * ENTER “1° FOR B SCANS  '2" FOR 4 SCANS® P 

650 9 77$=";INIFRNAL,4,1" | USES LEFT DISK DRIVE FOR RAW DATA FILES 
bh) CALL Readata 

470 CALL Shift 

680  fEEP 

670 Screen; | 

700 PRINT CHRE(12) ICLEARS SCREEN 

719 P1$="YOLTAGE® 

720 P2S="DIODF NUMER" | STRINGS FOR PLOTS 

730 PRS="YALTAGE vs, DIODE" 

740 P4S="RAW DATA" 

750 5=1024 

751 MAO, 

760 Xt=64 1 SET UP VALUES FOR PLOTS 

770 Yt=.04 

780 8 =Xm=4 

790 = -Ya=5 

791 = Xx=0 

B00 CALL Plot(Y2(#)) 1 DRAWS AXES ,ETC,.. 

810 GSTORE Grafts) | STORES GRAPHICS DISLAY JUST MADE 
820) CALL Dataplot(R,Y1(#) ,H) IPLOT NO-PARTICLES DATA 

839 PRINT * HO-PARTICLES PLOTTED ® 

840 PRINT "* 

Q50 PRINT * HOW MANY SCANS SEEM TO RE GOOD 2" 
B60 «= INPUT J 


71 


| 


Bae BR iNT WHICH SCANS ARE GOOD ?...1e...1,2,4,5,7, INCLUDE LAST CO 
MMA ILAST COMMA IS REQUIRED OR YOU HAVE TO HIT CONT.. TWICE 
B90 INPUT Scans(#) 

960 PRINT CHR$(12) 

918 PRINT USING "/////////" 

920 PRINT * AVERAGING THE SELECTED SCANS" 

938 GRAPHICS OFF 

940 = CALL Average(J ,Scans(#) , Y1i(#), Avi(#)). 

950 PRINT USING "@" 

960 GLOAD Graf(#) ILOADS GRAPHICS ARRAY RATHER THAN WASTE TIME RE-DRAWING 
970 GRAPHICS ON 


9780 CALL Result (Avt(#), X(#),H) ! PLOTS 1024 ELEMENT ARRAYS 
990 PRINT USING °//* 
1000 PRINT ° Average Intensity No~Particles® 


1010 ON KEY 0 LAREL "AVERAGEL® GOTA Screen 
1028 ON KEY 1 LABEL "PLOT-PARTICLES* GOTO 1060 


1030 PRINT ° PRESS KEY $ 8 TO RE-AVERACE OR # 1 TO CONTINUE" 
1040 Standby: 1 MANUAL CALLS THIS INTERRUPT DRIVEN PROGRAMMING 
1050 GOTO Standby 1 LOOPS, WAITING FOR USER TO DECIDE 
1060 PRINT USING "@* 

1070 OFF KEY 1 IHELPS AVOID CONFUSTON BY CLEARING THAT FOX 
1080 CALL Dataplot(B, Ye(#),H) 'PLOT PARTICLES DATA 
1096 PRINT USING °/* 

1108 PRINT ° PARTICLE DATA PLOTYED® 
1110 PRINT * FOR A HARD COPY OF THIS RAW DATA PRESS KEY $ 6" 

1120 PRINT * * 

1130 PRINT ° OR TO CONTINUE PRESS KEY # 1" 


1140 ON KEY 6 LAREL "HARD & RAW COTO Raw 
1150 ON KEY 1 LABEL "CONTINUE " COTO Select 
1180 GOTO Standby 


1170 Raw: CALL Plot(¥2(#),1) | THE ONE (1) IS AN OPTIONAL PARAMETER 
11890 CALL Dataplot(B, ¥2(#) ,H) | WHICH IS USED TO GET HARD COPIES 
1190 Select:! 

1200 OFF KEY 6 

1210 OFF KFY 1 


1220 PRINT CHRS(12) 
1230 PRINT USING *////* 


1240 PRINT * HOW MANY SCANS SFEM TO BE GOOD ?° 
1250 INPUT J 
1260 PRINT " WHICH SCANS ARE GOOD 7...1e,..1,2,4,5,7, INCLUDE 


LAST COMMA" 

1270 INPUT Srans(*) 

1280 PRINT CHR${(12) 

1290 PRINT USING "////////7° 

1300 GRAPHICS OFF 

1310 PRINT ° AVERAGING THE SELECTED SCANS* 
1320 CALL Average(] ,Scans(#), Y2(#) Av2(#)) 

1330 PRINT USING "8" 

1340 P4$="AVERACED SCANS* 


TZ 


2 & me eee 


1350 
1360 


1370 
1380 
139 
1400 
1410 
1420 
1430 


CALL Plot(Auv2(#)) 
GRAPHICS ON 


GSTGRE Graf(*) 

CALL Result(Av2(#) ,X(#),H) 

PRINT USING *//° 

PRINT " Average Scattered Intensity” 
ON KEY 2 LABFL "AVERACE2" COTO 1980 

ON KEY 3 LABEL "SURTRACT" GOTO 1450 

PRINT * HIT KEY #2 TQ RE-AVERACE OR KEY # 3 [0 


CONTINUF® 


1449 
1450 
1469 
1470 
1439 
1499 
1590 
1510 


GOTO Standhy 

FRINT USING *e" 

CFF KEY 3 

MAT Avi= Av2-Avt ISUBTRACTS NO-PARTICLES FROM PARTICLES 
GLOAD Craf(*) 

CALL Result(Avi(#) ,X4#) ,H) 

PRINT USING "///" 

PRINT " Plot of the Difference Ketween Particles and N 


o-Particles” 


1520 
1920 
1540 
190d 
1260 
1370 
1280 
1590 
1600 
1610 
182} 
1630 
1649 
1650 
1640 
1679 
1680 
1690 
1709 
1710 
1729 


1730 


PRINT * KEY #6 FOR BARD CePy’ 
PRINT * KEY #1 TO CONTINUE " 
ON KEY 6 LAREL ° HARD AVERAGE” GOTO 1570 

MN KEY {| LAREL " FILTER ° GOTO 1600 

COTO Standby 

CALL Pint(Avl(#), 1) 

FALL Resul t(Aui(#) ,X(#),H) 

PLOTTER 15 3, “INTERNAL” 

PRINT USING "//° 

PRINT * ENTER THE NUMBFR OF TIMES YOU WISH Tt) APPLY” 
Se THE DIGITAL FILTER FCR SMOOTHING" 

PRINT * EXAMPLE XHex KKH 10" 

PRINT " TAKES ABOUT 1.5 MINUTES® 

INFUT ° ZERO 1S ALSO ACCEPTARLE" Fil 

IF Fil=0 THEN Causs 

P4$="FILTERED DATA" 

CALL Filter(Ay1(#) Fil) 

PRINT CHRE( 12) 

CALL Plot(Aui(#)) 

TF (itd THEN MAT X= Ti 

CALL Result(Ayl(#) ,X(#),H) 

PRINT * Plot of the Difference Between Particles and N 


o-Particles" 


1749 
TER" 
iw 

1769 
1770 
1780 
1790 
1800 
1819 


PRINT ° AFTER APPLICATION OF A DIGITAL FIL 


PRINT * for a HARD COPY " 
PRINT * PREPARE the PLOTTER and" 
FRINT ° PRESS KEY € Q” 

PRINT * OR* 

PRINT " PRESS KEY # 9 TO NORMALIZE” 
ON KEY 8 LAPEL "HARD FILTERED® GOTO 1839 

ON KEY 9 LABEL " NORMALIZE * COTO Gauss 
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B20 
830 


1846 
1858 
1850 
1870 
1880 
1870 
1891 
1892 
1900 
1910 
1920 
1930 
1949 
1950 
1960 
1970 
1980 
1990 
et ng 
2010 
2920 
C030 
2031 
2040 
fall 
2060 
2070 
2080 
ON9§ 
e140 
126 
21 30 
P44g 
2150 
2160 
2170 
2180 


GOTO Standby 
CALL Plot(Avl!#), 1) ! The 1 15 for hard copy 
CALL Result (Avl(#) ,X(*) ,H) 
NFF KEY 9 
Gauss: OFF KEY 2 IKEY 1 OPTION IS STILI VALID 
OFF KEY 8 
ON KEY 1 LABEL ° FILTER " GOTO 1600 ! ALLOWS ONE TO SMOOTH OLD BATA 
1 WHEN REVIEWING REDUCED FILES 
GBid=1 ! THIS ALLOWS FILTERING EVEN AFTER ENTERING GRAPHICS MODE, THE 
1 DIODE COCRDINATES AREN’T NEEDED ANYMORE *SEE LINES NEAR CALL FILTER 
PRINT CH8$(12) 
'Av2 Array 15 Normalized BUT Avl(#) is Stil] Sayed For Re-werk if Needed 
MAT Ave= Avi =! TQ BEGIN WITH THE ARRAY IS ASSUMED TO 8E NORMALIZFD 
1 AND THE USER CHANGES THIS WITH INTERACTIVE GRAPHICS 
| IN THE SUBROUTINE “Compare” 
PiS="FILTERED VOLTAGE : 
P2S="THETA (rad) ISTRINGS FOR PLOTS 
P3$="INTENSTTY vs, THETA" 
P4$="CURYE FIT RESULTS" 
M5=,070 ISET UP VALUES 
Xt=.005 IFOR PLOTS 
Yt=.1 
Xn=2 
Yu=2 
Mew 
CALL PlattAyd(#}) 
MALL Result (Ayt(#), Ti(#),H) IPLOT OF NORMALIZED INTENSITY vs. THETA 
GSTORE Craf (x) | SAVES THE SCREEN IMACE 
CALL Compare(ful(#) HM MI Graf (*)) 
GRAPHICS (FF 
ON KEY 4 LABEL "OTHER APRAY" GOTO Choose 
ON KEY § LABEL ” TWO-ANCLE" GOTO Ruchele 
ON KEY & LABEL "STORE DATA" COTM Storedata 
OW KEY 7 LABEL * QUIT ®° GOTO Quit 
CALL Menvel ' PRINTS OPTIONS ON THE SCREEN 
GATO Standby 
Ruchele: |! 
CALL Twoangle 
CALL Menyel 
B9T9 Standby 
Storedata:! 
CALL Stare 
CALL Menve!} 
GOTO Standby 
Quit: | 
END 
SUB Averane(J ,Scans(#) ,¥(#) Av(#)) TAVERACES SFLECTED SCFNS 
MAT Ay= (9) (INITIALIZES THE ARRAY LOCAL TO THIS ROUTINE 
FOR J=1 TO J STEP j 
K=(Scans(1)-1)*1024+1 ! THIS COUNTER IS THE BEGINNING 
L=Scans(i)#1024 1 ANT THIS ONE THE END OF BLOLKS 
FOR Ti=K TOL 1 OF 1024 IN THE OVERALL DATA 
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L2={{-K+1 1 THIS COUNTER IS ALWAYS DETWEEN 
Av(12)=Ay(12)+¥ (11) \ Mp ite 


2539 { 

2740 NEXT II 

2350 NEXT I 

2550 MAT Av= Au/(J) 

2370 SUBEND 

2580 SUR Dataolot(R, ¥(*),H) 'R IS 4 OR Q(THE NUMBER OF SCANS) 
2390 COM /Max/ M7,M5,Xt,¥t,Xm, Ya, Xx 

2400 LDIR 0 

C410 = LORG 4 

2420 FOR J=1 TN RB 'EACH SCAN 

2439 MOVE S00+J*50,47-,05 | SEE NEXT LINE 

C449 LAPEL J; | HELPS KEEP TRACK OF FACH SCAN AS IT APPEARS 
2450 LINE TYPE 1 ’ 

2440 K=(J-1)*102441 IRECINYING OF EACH SCAN ARN 

2479 L=J¥{n24 ITHE END WETHIN THE TOTAL BLORK 
2480 MOVE 1,Y{K) 1 MOVE TO THE FIRST PAINT 

2499 FOR (=K TO L STEP H 

2000 X=1-K+] 1 THIS GIVES 1 TO £924 FOR ARSCISSA 
cold PLOT X,¥¢1) 

eved NEXT I 

253) NEXT J 

2340 PENLIP 

eof MB=4) 

cu SUPEND 

25b8 SUP Resul t(Tav(#) ,D(*) ,H) ITHE AVERAGE INTENSITY IS PLOTTED 
2u70 CLIP ON 

2588 MOVE DC), Tay) 

co90 FOR I=1 TO 1824 STEP H 

2HO0 DRAW DIT) ,Jav¢1) 

CA1O NEXT I 

cb20 PENUP 

2430 SUBEND 

eb3t Snit: |! 

C640 SUB Plot(¥(*) ,QPTIONAl. Hard) ITF THE OPTIONAL (Hard)IS PECEIVED 
2650 COM /Max/ M7,MS,Xt, Vt, Xa, Ya, Xx ITHE FIGURE GOES 10 THE PLOTTER 
2660 COM /Plots/ P1$(20) ,P2sl20) ,P3$(20) ,P4$(20) 

C70 MP=MAX(Y(#)) #1, 1 1 A SCALING VARTABLE 

2680 Mé=,1*H7 | ANOTHER SCALING YARTARLE 

e700 GINIT IGINTT IS JUST C9OD PRACTICE SU YOU KNAW 
e710 IWHERE YOU ARE RECINNING 

C720 «SELECT NPAR ITHIS DFIECTS IF FHE HARD COPY TS DESIRED 
2730 CASE 1 'IT COULD BE DONE WITH IF THEN LOGIC 

2740 PLOTTER IS 3," INTERNAL" IRUT 15 PRESENTED FOR FAMILIARI ZATION 
2798 CASE 2 ISINCE IT IS HORE POWERFUL. IN COMPLEX 
2740 PLOTTER IS 705, “HPGL" ISITUATIONS 

2770 END SELECT 

2730 GRAPHICS ON 1 TO BE ABLE TO SEE THE PLOT 
eyu CSlit 9,.6 

2600 DEG 

CB10 LDIR 0 
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erm 5 wanes ee ee 


onan - ee ee me 


ead 


2920 LORG 5 

2830 MOVE 75,95 

2840 LABEL USING "K";P4$ 

e850 MOVE 75,90 

2860 LABEL P3$ 

2870 MOVE 75, 20 

2880 LABEL P2$ 

2890 LDIR 90 

2790 MOVE 18, 60 

2910 LABEL P1$ 

2720 VIEWPORT 30,125,390, 1 SCREEN UNITS FOR MARGINS 
2930 FRAME | DRAWS A BOX 
2940 WINDOW Xx, HS, -MS, 4? 

2950 AYES Xt,M6,Xx,0,Xa,1,9 

2960 LDIR 0 

EU VUTSIZE 33a 

2980 LORG 8 

2590 CLIP OFF 

3000 FOR J=-M4 10 M7 STEP M7/10 INUMRER THF Y AXIS 
3049 MOVE Xx ,T 

3020 LAREL USING "#,MD, DDDD";1 

3039 NEXT I 

S00 Vale 3,06 

S050 LDIR 70 

3069 LORS 8 

3070 FOR J=Xy TO MS STEP Xt INUMBER THE X AXIS 
3080 MOVE 1,-Mé 

5090 LAREL USING "#,K";J 

3100 NEXT | 

5410 PENUP 

3120 SUBEND 

3130 SUB Compare(Avi(*),H MME INTEGER Graf (*)) 
3149 COM /Gauss/ TtCe),L 

3150 COM /Hrdgaus/ Av2(#) 

3168 COM /Max/ M7,MS,Xt,¥t, Xa,Ya, Xx 

3201 Poot: | 

S202 PRINT " INPUT REFERENCE ANGLE " 

3293 INPUT Tt 

S205 FOR N= 10 1024 STEP 1 

3296 IF TI(N))=Tt THEN 3208 

3207 NEXT N 

3208 Tt=f 

3299 FOR J=N-5 TO NtS 

3c 10 It=It+tAvl (1) 

Self NEXT J 

Seve Al=It/11 


3220 ORR KeEY IGETS RID OF ALL LABELS ON KEYS 
3233 ON KEY 3 LABFL "MENUE™ GOTO 4169 
3240 PLOTTER IS 3,"INTERNAL" | IN CASE A HARD COPY WAS JUST MADE 
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3290 
3250 
5270 


Hard=0 1 (0) $0 ONE DOESN] EX 
pameeline=! TE THT TEAL HORMELTS TW VAL Lie 
D=19, | INITIAL PATICLE MEAN DIAMETER IN MICRONS 


3280 Change: CSIZE 4,.6 


3290 
3300 
3310 
3320 
Joal 
3349 
350 
3379 
oono 
3780 
3440 
3430 
3469 
347 
SIE” 
3480 
3458 
5509 
$510 
Rhy!) 
3539 
a54n 
3541 
$558 
3569 
357) 


B80 Wait: 


3598 Puls 
36.6 
3510 
3625 
3530 
3648 
3699 
Sbb8 
3679 
3680 
3699 


PRINT USING "//////* 

PRINT * STAND BY FOR CURYE" 
MVE 0,1 

CLIP ON 

DIM T(1:1824) 

FOR I=1 10 1024 STEP H 

TEL) 201 EXE (-D*24C (TH CT) 42-1142) eC (M1, S7HPT)/L)42) )) 

DREW THT), TCT) 

IF (TCL) *( (PIRD#NE /L)) 93,0 THEN 3440 


NEXT | 

PENIIP 
TF Hard THEN 355 1EXITS ROUTINE IF A HARD COPY WAS JUST HADE 
PRINT USING "///////1/" 
PRINT * USE THE KNDR TO VARY THE PARTICLE 
PRINT ° QR HIT KEY # & FOR HARD COPY" 
PRINT * " 
PRINT * KEY #9 ALLOWS YOU 1" 
PRINT * SELECT NEW REF ANGLE" 
PRINT > * 
PRINT * HIT KEY #3 70 GET OUT" 
POINT USING “//////" 
QFF KEY 0 


ON KEY 6 LAREL "HARD COPY" COTO Hardgavss 
ON KEY 9 LABEL “NGRMALIZE” GOTO Poot 


QN KNOPR .05 COTO Pulse 1(,65) IS TIME INTERVAL IN WHICH 
GOIN Wart 1 PULSES FROM THE KNOW APE 

e: PRINT WSING "e° 1 COUNTED AND THIS NUMPER IS 
Strani$="Mean Diameter " 1 USED BY THE INTERRACTIVE 
Strnq2$=" Microns" 1 GRAPHICS TH VARY THE PARTICLE 

Count=KNOBX 1 SI7E AND PLOY THE ASSOCIATED 

D=DROUND(D+Count/15, 2) 1 GANSSTAN APPROXIMATIONS OF 

GLOAD Graf(*) 1 SCATTERING PROFILES 

LDIR 0 

LORG 8 


HOVE KS, .8*M7 
LABEL Strngi$SAVAL$(D) &Strnges 
GOTO Chance 


3700 Hardgauss: ! 


3719 
3720 
3739 
3740 
3730 
3768 


PRINT " PREPARE the PLOTTER" 
PRINT * PRESS CONTINUE for" 
PRINT ° HARD COPY" 

PAUSE 

Hard=1 1 SO THAT SUREXIT OCCURS AFTER HARD COPY 
CALL Plot(Av2(*),1) 1 TEXHARD COPY 
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377 
3780 
3790 
33809 
3016 
3820 
3849 
4160 
4170 
4180 
4190 
4206 
4210 
4220 
4230 
424) 
4250 
4260 
A279 
4280 
429) 
4300 
A310 
4320 
43g 
4340 
4350 
4360 
4378 
4390 
4370 
4400 
44}) 
4429 


CALL Result(Ave(), TH(#),H) 


PRINT USING °@° 
LORG 8 
MOVE M5, .8*M7 
LDIR 0 
LABEL Strrg1$AVAL$(D)A&Strng2s 
COTO Change 
Subexit: 
SUBEND 
SUB Menvel 
PRINT USING "@" 
PRINT " YOU CAN RE-NORMALIZE" 
PRUNT GOGH GUOBUEBOB GU HOHE GUGUE OH OOOO BHBOUGUD OG CY! 
PRINT * YOU CAN RE-AVERACE (New Data finiy)" 
PRINT “GHGRH HEB BEBE HO EURO COG OGEHOOURORO OEY 2220008 bE ZEX x" 
PRINT ° PRESS KEY # 4 TO LOOK AT OTHER DATA® 
FRINT ° OLD OR WFW 
PRINT ° MOTOR / EXHAUST® 
PRN T° 3636 38 96 96 06 360636 963098 38 309 9 EE 2 SME SE SEE IE CEI A 
PRINT * PRESSIKE) Seams 
FRINT ° FOR THE ’ TWO-ANGILE ’ METHOD” 
PRINT ° OF PARTICLE SIZING” 
PRURMLESISIELIC ETC ETT P TT LISI ST CELE Li OStl Sct ea titsetet eel ti ett © a 
PRINT " TO STORE THE REDUCED DATA PRESS KEY ¢ & F 
PRINT UGQHGHBHBUBOHE BH OBO URE OUR UDO GOBEXUERO HEE HEHES x # 
PRINT * TO END THIS SESSION PRESS KEY # 7" 
SUPEND ; 
SUB Displayl (Old) 
PRINT ° TO REDUCE NEW DATA PRESS KEY # 1° 
PRINT °° 
PRINT ° TO REVIEW PREVIOUSLY REDUCED DATA PRESS KEY # 2" 
PRINT " ° 
QN KEY 1 LABEL ° NFW DATA’ COTO New 
ON KEY 2 LASEL " OLD DATA" GOTO Review 


4439 Wait: COTO Wait 


44a 
44%) 
446) 
4470 
4480 
449) 

290 
4310 

20 
4530 
454g 


New: PRINT CHR$(12) 

PRINT ° PUT THE DISK IN THE LEFT DRIVE IF IT IS NOT ALREANY® 
PRINT °* 

PRINT ° ENTER THE NAMES OF THE FIRST AND SECOND FILES." 
PRINT "° 

PRINT ° EACH FILE HAS DATA FROM BOTH DIODE ARRAYS," 
POINT "* 

PRINT ° FIRST 1S NO-PARTICLES & NEXT IS PARTICLES ----Di$,D2$° 
PRINT "2 

Old=0 

SUBEX IT 


4350 Review: Qld=1 ITHIS VARTARLE IS PASSED T!) THE MAIN PROGRAM TO INDICATE 


4°60 
4570 


'THAT THE DATA TO RE READ TS ONE (1) SCAN AND THAT 
ITHE AVERAGING ROUTINES WILL NOT BE USED 
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4°30 
490 
4600 
4510 
462) 
4630 
464) 
4650 
4669 
4670 
4680 
$699) 
4700 
4710 
4720 
4739 
4740 
$700 
4760 
4779 
4780 
4790 
4610 
4310 
4820 
4830 
4840 
4850 
4560 
487) 
4890 
4399 
4900 
4919 
4920 
4930 
4974) 
A950 
496) 
4979 
4980 
4790 
1000 
3O10 
2028 
1039 
3040 
30590 
30.60 


SUPEND 

SUB Readata 
COM /Readata/ B,P,H,Q3${20) 48120) 22920), Yt) BUFFER ,Y2(*) BUFFFR 
Xyz=] FILE POINTER VARTARIE 
#=8 | NUMBER OF SCANS IN MOTOR DATA 
R1=65534 INUMBER OF BYTES MF MOTOR DATA 
IF P=2 THEN Xyz=4097 IRECCRD $ WHERE 4 SCANS DATA RESINS 
IF P=2 TREN R1=32768 INUMRER OF BYTES OF 4 SCAN DATA 
TF P=2 THEN R=4 INUMPER OF SCANS 
PRINT? = * 
PRINT ° READING DATA FROM FILE ON DISK’ 

ASSIGN @Diskt TO Q38&2Z2% 

ASSIGN @Pisk2 TO Q4$8772$ 1OPEN 1/0 PATHS 


ASSIGN @Buffl TO BUFFER Y1(*) 
ASSIGN #Ruff2 TO BUFFER Y2¢#) 
CONTROL @Diskt ,osXyz 


TRANSFER @Disk] 10 @Buff1;COUNT R1 IRFADS NO-PARTICIE DATA 
WAIT FOR EQT @Diskt | NAT AN OVERLAPPING TRANSFER 
CONTROL @Disk2,o;Xyz IXyz SETS DISK FILE POINTER 10 

VEITHER MOTOR OR EXHAUST DATA 
TRANSFER @Oisk2 TQ @Buff2;COUNT Ri 'READS PARTICLE DATA 


WATT FOR FOT @Diske 
ASSISN @Diski TO * 


ASSIGN PDiske TO * HUST GOOD PRACTICE 19 CLOSE 
ASSICN @Ruffl TO % 11/0 PATHS 
ASSIGN @Buff2 TO * 

SUBEND 


SUB Plot2(F D(x) ,X) 


DIM Titleit${25],Theta$(20) 

GINIT 

IF X=1 THEN PLOTTER 1S 705, *HPGL" 

DE IDEGREES for LAREL DIRECTIUIN 

ERAPHICS QN 

VIEWPORT 35,125,35,85 

Max=10#INT( (MAX (D(C) )410)/10) 

WINDOW 0¢1,2),D(E,2) ,0 Max 

G=(D(E,2)-DUL ,2))/(E-1) #4 1AN X GRID LINE EVERY 4th POINT 

F=(INTCE/4)-1) #4 

IF F=0 THEN F=4 

CLIP D(4,2)-G,D(F ,2)+286,0 Max IMAKES GRID UNIFORH 

GRID 6,5,0(4,2)-6, 8 

CLIP OFF 

LCRG 8 

LDIR 90 

CoE 4.3 

FOR 1=4 TO E STEP 4 IPUTS NUMBERS ON X AXIS 
MOVE D(T,2),0 
LABEL USING ".DDDD";D(I,2) 

NEXT I 


Ts, 


1488 
0399 
5400 
v4!) 
9429 
7439 
5440 
u450 
7460 
v470 
3480 
3499 


Lhe 

FOR I=10 10 Hax STEP 10 
MOVE D(4,2)-6,,1 
LABEL USING "#,K"s1 

NEXT I 

CSITZE 6,.6 

Titles="THO-ANGLE METHOD” 


Titlei$="For Various Angle Ratios” 


Sizes="SiZE (microns)" 
Thtag="THETA (rad)" 
Sub$="1" 

LOIR 90 

LORG 5 
R=D(4,2)-C-(D(E,2)-DUI ,2))/10 
MOVE KR Max/2 

LAREL Size$ 

LDIF 8 

A=(D(E, 29401 ,2))/2 
MOVE A,-Nax/4 

LABEL Thta$ 

LORG 3 

MOVE A,-Max/4 

LABEL Seb$ 

LIC 5 

MOVE A Maxal, £5 
LAREL Titles 
(Sea 

MOVE A Max#i.05 
LABEL Titlet$ 

PENUP 


SUEEND 
SUB Distribution(D(*), Tratio,£) 


LINE TYPE 1 
CLIP CN 
Tratiot=VAL$(Tratio) 
MOVE O(1,2),00,1) 
FOR I=1 IDE 

DRAW D(T,2) DUI, 1) 
NEXT I 
LORG 4 
LINE TYPE 1 
COIZEM An. 3 
HOVE D1 ec) biel lite 
LABEL Tratio$ 
PENUP 


5058 SIREN) 
JJ49 SUB Twoangle 


INUMPER Y AXIS 


ISTRINGS FOR PLOTS 


| ANGLE RATIO 


ITHIS SUBROUTINE PLOTS THF 
IPARTICLE SIZES DERIVED USING 
IVARTOUS ANGLE RATION APPLIED 
'TO THE DATA OVER A RANGE OF 
TANGLES 
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Jou0 


Peete Ce LAS SCS ete lire etree sites eee e tert ses cece sete eer artseses: 


De LiMBRTSSISC SST LS 22 SUBPROGRAM ‘’TBUANCLE’ KK EERE EEE ERE EEK YX 
9570 | RRRRRRRE EERE ER PARTICLE SIZING BY eet e settee eases eS: 
Wil MES eS ctSesetees 3: MEASURING THE RATIO XCEL LEE EE EE LEX 
WWMM CECLOLSS tetas ¢ | QF INTENSITY AT Pe eSetiLetTteett ele: 
THON | EEN ERHEX EXE REE TWO ANCLES HK ELRE NY EERE CREE EEK 8 
WIDER LISI SEC CLC s sts c ets e tite t ects t tse stat eseci ite tess tt tete tects 


9620 OPTION RASE { 


5630 


9450 


5660 DIM D(200,2) 


y570 PRINT 


COM /Two/ Avitife+) ,M,Al,F 
340 COM /Gauss/ T141024),L 


ITHE GAUSS IS NOT USED HERE. 
VEY THE COW RLOCK HAS THETA 
|AND WAVELENGTH 

CHR$(12) 


S680 PRINT USING "///////" 


JE90 PRINT 
1700 PRINT 
O710 PRINT 
u7ed PRINT 
0730 PRINT 
3740 PRINT 
S750 PRINT 
2760 PRINT 
3770 PRINT 
780 PRINT 
9799 PRINT 
2900 PRINT 
J810 PRINT 
otal PRINT 
3830 PRINT 
W840 INPUT 
850 
wef} 


OFF 
X= 
2870 Regin: 


. THIS SURPROGRAM USES, THE TWO-ANGLE HETHOD DESCRIBED RY BUCHFIE* 


. TO CALCULATE PATICLE SIZE FOR VARTOUS ANGLE RATIOS AND ANGLES, ° 


IT 15 HOPED THAT THE CURVES WHICH RESULT WILL SHOW ° 


WHICH SIZE IS THE MOST PROBABLE. ° 


Li) 


"AFTER NOTING FROM THE RAW DATA WHICH ANGLES CONTAIN THE CFNTER LOBF" 
: ENTER Re RReMRRRERERTHE SHALLEST USEABLE ANGLE," 
: THE SMALLEST ANGLE RATIO,AND" 
THE STEP SIZE RETWEEN ANGLE RATOS* 
YOU WISH TO EXPLORE’ 
MWilepive,. 4 
thetal, theta-ratio, theta-step’,Q,A,B 


3 EXAMPLE 4% HH HER EEE ER 

" ENTER 
KEY 

IGRAPH WILL APPEAR ON SCREEN 

IRUNNING CONTINUES HERE WHEN A HARD COPY JS DESIRED 


JEB0 Wan De Hulst and Gumprech & Sleepevich exolain that the chanoe in 
1890 twavelenoth of the beam is accounted for by dividing by the index 

3919 lof refraction of the medium, €#K#XEXXL=L/M] EXER KEE 

J710 C=(L/MI/.57/PT) “2 | a see page 15 of nasa terh paoer 215 


5928 FOR N=1 TO 1024 

IF Q¢TICN) THEN 5990 

2940 NEXT N 

9950 FOR Tratio=A TO 3 STEP R 

FOR J=N TO 1010/Tratio STEP 10 


53930 


5950 
3970 
598e 
3790 
6009 
5019 
6020 


It9 see this 15 a convenient constant 
IFINDS POSITION OF MINIMUN AXLE 


IWARTOUS ANSLE RATIOS 
ISCETS THE RANGE OF POSSIBI_E 
| ANGLES (THETA 1) 
Thtal=fity) 
Thtad=Th(J*Tratio) 
Deltatheta=Ihta2*e-Thtal “2 
11=0 
12=0 


! THETA 2 
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Mh RRA, Rd PSN Yi 


6050 NEXT I 'OF A STEADIER CALCULATION 
6060 T1=]t/ti 

6070 FOR I=INT(J¥Tratio-5) TO INT(J#Tratio+5) 

6080 c=l2tAvi (1) | DOME HERE FOR THETA2 FOR 
6090 NEXT I 1 THE GIVEN ANGLE RATIO 

5100 l2=T2/i 

6110 E=(J-N)/10+1 ITHIS IS A COUNTER FOR THE ARRAY CONTAINING 
6120 'PARTICIE SIZE AND THETA! FUR THE GIVEN 
6130 IANGLE RATIO 

6140 IF TidI2 OR Ti=0 OR T2t=0 THEN DIE, 2)=8 ALLOWANCE FOR IF THE DATA 
6159 IF [idI2 OR Tis0 OR [2¢=0 THEN GOTO Xcomp '1S NOT WELL BEHAVED 
6160 D(E,1)=SQR(-C/Deltatheta eLOG(1°/11)) (TWO ANGLE METHOD 
6179 'VALUE OF DIAMETER based on INTENSITY RATIO 
6189 [For a given ANGLE RATIO and ANGLE THETA! 
6190 Xcomp: D(E,2)=Thtal 

6200 Spar=PI¥D(E, 1) *M1/L ITHIS IS PI¥D/LAMRDA -- THE SIZE PARAMETER 
h10 Thar=Snar#thtae ITHETA PAR FAR THE LARGER ANGLE 
6220 ~ TF Thar)4 THEN J=1010/Tratio (THIS ENDS JHE DO LOOP FAR THIS 
6230 IF E=1 AND Thard)3 THEN 6330 IANGLE RATION SINCE THE GAUSSIAN 
6240 '!the above line ends all calculation if 11S NOT VALID WHEN Thar ) 3 


6250 ‘the first element (E=1)failed the test 

62769 NEXT J 

6270 | FIRST TIME THRQUGH-- D HAS THE MOST ELEMENTS IT WILL HAVE 
628 9 | =S0 SET UP THE GRAPH USING D’s PRESENT PARAMETERS 
6290 PRINT CHR$(12) 

6300 IF Tratio=A THM CALL Plot2(E ,D(#) ,X) 

bat 8 CALL Distribution(D(#) , Tratio,e) 

6320 NEXT Tratio 

6330 ON KEY 2 LABEL "HARD COPY” COTO Hard 

6240 ON KEY 3 LABEL “MENUE " GOTO Subexit 

6350 PRINT USING °/////" 


6360 PRINT ° YOU CAN GET A HARD COPY BY PRESSING KFY # 2° 
6370 PRINT " OR EXIT THIS ROUTINE BY PRESSING KEY ¢ 3" 
6380 Standby:GOTO Standhy 

6390 Hard: X=1 '§ VARIABLE TQ CONTROL PLOTTER 


6400 COG Ceqin 

6410 Subexit: GINIT 

642) GRAPHICS OFF 

6430 SUBEND 

6440 SUB Shift 

6450 COM /Readata/ B,P,H,Q3$f20) ,Q4$(20) 2278120) ,Yi(#) BUFFER, Y2(#) BUFFER 
6460 DIM E(1:8192) 

6470 PRINT CHRS(12) 

6480 PRINT USING "//////7" 

6498 PRINT " RAW DATA IS BEING SHIFTED TO CORRECT FOR SMALL GAPS BETWEEN SCANS* 
moO) PRINT °* 

6510 PRINT " BE PATIENT, IT’S A LONG SET OF DO-LOOPS’ 
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6020 
6539 
6548 
b550 
6568 
6570 
658) 
6590 
6600 
6619 
6420 
6538 
6640 
6450 
6660 
6670 
6680 
5659 
6709 
6710 
6720 
6730 
6749 
6750 
4769 
677) 
6780 
6799 
6800 
6810 
6220 
6830 
634) 
6850 
6869 
5870 
6880 
6870 
6909 
6910 
692) 
6930 
6949 
6950 
6960 
597) 
6980 
6590 
7008 


| THERE ARE SOME SMALL GAPS a CANS AN 0 Bhar re want 


| SHIFTS THE DATA SO THAT THE FIRST DIODE DA 
| TO THE VERY BEGINNING OF ITS 1024 BLOCK IN THE OVERALL ARRAY 
THE FIRST SET IS RIGHT ON, THE NEXT 159 ONE OFF,THE THIRD [5 
| THO OFF ,,50 FORTH. THIS MAY NOT BATTER WITH OUR RESALUTION 
AND 1S PROBABLY DUE TO THE MEMORY CARD CYCLING AT THE END OF 
| SOME SCANS, SEE THE CIRCUIT TIMING DIAGRAM IN THE THESIS. 
SELECT B 
CASE 4 'EXHAUST DATA HAS 4 SCANS 
K=3 
CASE 8 IMOTOR DATA HAS 8 SCANS 
M=7 
EXD SELECT 
FOR K=1 10 2 
IF K=1 THEN HAT E= YI IOPPERATES ON NO-PARTICLE AND PARTICLE SETS 
IF K=2 THEN NAT E= Ye 
FUR J=0 10 4 
IF M=7 AND J(=3 THEN 6800 'One 4096 Block Doesn’t Heed Shifting 
FOR T=(J)#1024+1 TO (J+h)*1024 IRLOCKS OF 1024 
IF 4=3 THEN L=I+J lArray B Has the worst 
Problem With Shifting Data 
IF M=7 THEN L=I+] lArray Dis Allways off by one 
IF Ld(J+1)41024 THEN L=(J+1 #1024 
1JUST TO AVOID PROGRAMMING ERROR AT THE END OF THE APRAY 
E(I)=E(L) ITHIS SHIFTS THE RATA 
IDEPENDING ON ‘L’ ARRIVED AT 
'BY MESERVING RAW DATA 
NEXT J 
NEXT J 
IF K=1 THEN MAT Yi= E 
IF K=2 THEN HAT Ye= E 
NEXT x 
CUBREND 
SUR Store | THE LARGE ARRAY CONTAINING MULTIPLE SCANS HAS BREEN REDUCED 
| 70 A MEAN SCATTERING PROFILE BY AVERAGING AND FILTFRING, IF 
| YAU FEEL THAT STORING THIS DATA I5 NECESSARY THIS RMUTINE 
1 DOES 17, IT SAVES DISK SPACE TO STORE THE RESUIIS THEN 
| ELIMINATE THE RAW DATA IF YOU FEEL CONFIDENT THAT THE 
| REDUCTION IS THE BEST THAT CAN FE. [N QTHER WORDS, 
| DA NOT PURGE A RAW DATA FILE UNLESS YOU ARE ABSOLUTELY SURE 
| YOI) WON’T WANT JO REDUCE IT AGAIN. 
COM /Hrdgaus/ Av2(*) lAv2(#) IS THE REDUCED DATA 
DIM Data(1:1024) SUFFER IBUFFER USED TO TRANSFER 10 DISK 
MAT Data= Ave 


PRINT CHRE(12) 

PRINT USING °//////" 

PRINT " A SUCCESTED METHOD OF NAMING REDUCED DATA FILES IS AS FOLLOWS" 
See. oN-===——- MOTOR BEAM ® 

SNS Sera EXHAUST BEAM" 
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7010 
7020 
7030 
7040 
7050 
7060 
7070 
7080 
7090 
7109 
7119 
7120 
7130 
7140 
7150 
7160 
7170 
7199 
7198 
7200 
7210 
722 

7230 
7240 
7250 
7260 
729 

7280 
7290 
7300 
7310 
7320 
7330 
7340 
7350 
7360 
7379 
7328 
7390 
7400 
7410 
7420 
7430 
7440 
7459 
7460 
7470 
7486 
7490 


PRINT © [------- CALIBRATION , IF NO ‘C’ THEN AN ACTUAL FIRING IS ASSUMFD" 
PRINT " G-------GLASS BEADS 

PRINT ° IF NN ’C’ THEN ’G’ STANDS FOR ‘GAP’ PROPELLANT® 
PRINT °  A-=-===-ALUMINUM OXIDE" 

PRINT * PPP-=-=-CHAMRER PRESSURE FOR RUN OR OTHER..,,(NOZZLE TYPE)" 
PRINT ° PARTICLE SIZE FOR CALIRRATION 

FRINT * MMDD----MONTH, DAY OF RUN OR CALIBRATION" 

PRINT ** 

PRINT * EXAMPLE:  HCAL25JN12" 

PRINT ** 


PRINT * MOTOR BEAM CALIBRATION USING ALUMINUM OXIDE FROM 1 TO 25 MICRONS" 
FRINT © QN JUNE 12" 

PRINT 
FRINT 
PRINT * EXHAUST REAM DATA FOR GAP PROPELLANT AT 350 psi ON AUICUST 10" 
PRINT PLACE A DISK IN THE RIGHT HAND DRTVE" 

PRINT " ENTER THE FILENAME YOU WISH TO USE FOR THIS REDUCED DATA" 
INPUT A$ 

CREATE BDAT A$,512,16 

ASSIGN @Disk TO A$ 

ASSIGN @Ruff TO BUFFER Data(*) 

CONTROL @Buff,3;1,8192, 1 IRUFFER IS FULL 
TRANSFER Shuff TO @Disk;COUNT B192 

WAIT FOR EOT @Disk 

ASSIGN PEuff TO * 

ASSIGN @Disk TO # 

SUBEND 

SUR Review(Ay1(#)) 

DIM Data(t:1024) FUFFER 

PRINT " EACH REDUCED FILE CONTAINS ONE SCAN OMLY. YOU HUST KNOW IF IT" 
PRINT " IS EXHAUST , MOTOR CAVITY , OR CALIBRATION DATA. * 

PRINT ** 


PRINT * THE DISK WITH THE REDUCED FILE SHOULD KE IN THE RIGHT-HAND DRIVE" 
PRINT °° 


EXAMPLE: XGSSOAUIO" 


PRINT ° ENTER THE FILENAME OF THE REDUCED DATA" 
PPENT ° TO BE REVIEWED" 

PRINT °* 

TNFUT AS 


ASSICN BDisk TO AS 

ASSIGN @Buff 10 BUFFER Data(*) 

CONTROL SDisk, 5st 

CONTROL @buff,3;1,0,1 ITHIS IS AN EPTY BUFFER 
TRANSFER @Disk TO @RuffsCOUNT 8192 11024*8=NUMSER OF BYTES 
WAIT FOR ENT BDisk 

ASSIGN @Disk TO x 

ASSIGN @Buff TO * 

MAT Avi= Data 

SUREND 

SUB Filter (A(*) Fil) 
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fo DIM E(1:1024) 
nt PRIN hibeet 
7520 PRINT USING *////////" 


7530 PRINT " FILTERING THE SCATTERING PROFILE" 

7340 FOR J=1 10 Fil 

7950 FOR [=1 TQ 1014 

7966 B=I+1 ITHIS 1S A Symetric Moving Average 

7570 C=I+2 ITYPE OF DIGITAL FILTER. EACH 

758) D=1+3 IDATA POINT IS EQUALLY WEIGHTED 
=90 F=[+4 'IN THIS CASE BUT THIS CAN RE 

7600 G=1+5 ICHANGED IF ONE DETERMINES THAT 
7610 H=1+6 IFEWER POTNTS WITH UNEQUAL WEIGHTS 

7629 K=[+7 WOULD BE FASTER OF GIVE PETTER 
7630 L=1+8 IRESULTS, THIS TYPE OF FILTER WAS 

7649 N=1+9 HISED SINCE IT INTRODUCFS HO PHASE 
7650 P=[+10 'LAG (Angular Resolution) . 


7660 ECG )=C ACT) FAC BR) 4A(C)4A(D) FACE) +A(G) HACH) FAK) FA(L) ACN) FAP) )/11 
7670 NEXT I 

7680 HAT A= F 

7599 NEXT J 

7700 = SUREND 
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